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Aims. The effect of frequency-dependent AGN core positions ("core-shifts") on radio Very Long Baseline Interferometry (VLBI) 
global astrometry measurements is investigated. 

Methods. The basic equations relating to VLBI astrometry are reviewed, including the effects of source structure. A power-law 
representation of core-shifts, based on both observations and theoretical considerations of jet conditions, is incorporated. 
Results. It is shown that, in the presence of core-shifts, phase and group-delay astrometry measurements yield different positions. For 
a core displacement from the jet base parametrized by Ax(A) = kA 13 , group delays measure a "reduced" core-shift of (1 -yS)Ax(A). 
For the astrophysically-significant case of /? = 1, group delays measure no shift at all, giving the position of the jet base. At 8.4 GHz 
an estimated typical offset between phase and group-delay positions of ~170yuas is smaller than the current ~ 250 //as precision of 
group-delay positions of the sources used to define the ICRF; however, this effect must be taken into account for future measurements 
planned with improved accuracy when comparing with optical positions of AGN to be obtained with the GAIA mission. 
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1. Introduction 

In astronomy the highest angular resolution achieved on a rou- 
tine basis is provided by the technique of radio Very Long 
Baseline Interferometry (VLBI). Observations with interconti- 
nental baselines can be used to image sources with sub-mas res- 
olution, and astrometry with a precision ~2 orders of magnitude 
smaller is possible. VLBI positions of a few hundred point-like 
extragalactic radio sources are used to define the International 
Celestial Reference Frame (ICRF; see Ma et al. 1998; Fey et 
al. 2004). Most radio sources are not points on this angular 
scale, however, and are far from circularly-symmetric, the emis- 
sion rather arising from collimated jets of relativistic plasma. 
Furthermore, different regions of the emission have different 
spectral indices. The question as to what points these positions 
refer to is not just of academic interest. The ESA space mis- 
sion GAIA (e.g. Lindgren et al. 2007), due for launch in 2011, 
will measure the positions of many AGN at optical wavelengths 
with precisions as good as 24 //as. Furthermore, proposed im- 
provements in VLBI instrumentation (Petrachenko 2006) hold 
promise for much improved VLBI global astrometric precision. 
Alignment of a GAIA-based optical reference frame with the 
ICRF will require a physical understanding of the location of 
the radio positions at this level. 



1.1. Astrometry using interferometer phase 

The high resolution of VLBI results from the fact that the visibil- 
ity phase, <p, is sensitive to changes in the geometric path length 
difference in the two arms of an interferometer. For a baseline of 
length L this difference is LcosO, where is the (instantaneous) 
angle between the source position, x, and the baseline direction 
in the sky. For a small position offset in the sky, Ax, the com- 
ponent in the resolution direction (the direction of arc on the 



sky), Ax , gives rise to a change in path difference of L sinO Ax 
resulting in an "astrometric" phase change 

A<p astr — (2nv I c)L sinO Ax 
where L sinO is the baseline length projected in the source di- 
rection, v is the observing frequency and c is the speed of light. 
The resolution is characterised by the lobe (or fringe)-spacing, 
Qiobe = c/ (vL sinO), the position change in the resolution direc- 
tion which changes <p by one turn. For baseline lengths of several 
thousands of kilometers and observing frequencies above 5 GHz, 
Sliobe can be less than 1 mas. Radio source images with resolu- 
tions (beam FWHM) ~ 0.4 Q/ /, e are made using VLBI arrays, 
and the separation between unresolved features in such images 
can be determined to a small fraction of the beamwidth. The 
relative position of two point sources close together on the sky 
can also be directly measured using the difference of their inter- 
ferometer phases. For measurements with high signal-to-noise 
ratios the precision can be a small fraction of Cliobe, ~ a few //as. 

However, <p is only measured modulo 2n whereas the ge- 
ometric path length difference, LcosO, can be hundreds of mil- 
lions of wavelengths. Furthermore, even with an accurate knowl- 
edge of instrumental phase and clock offsets, telescope positions 
and the orientation of the baseline in space, the signal propa- 
gation through the neutral troposphere and the ionosphere above 
the telescopes adds an unknown number of turns of phase, which 
varies on a short timescale. The additional path, A, through the 
troposphere (~2m at the zenith) is non-dispersive, adding an ad- 
ditional interferometer phase <p trop = (2nv/c)(A B -A A ) where the 
superscripts A, B refer to the two telescopes. However, the path 
through the ionosphere, /(v), is dispersive, producing a phase 
advance at each telescope, and an interferometer phase change 

4> ion = -(2tt/c)(A b - h A )v- 1 
where I\ is the path at unit frequency. For close source pairs 
(whose atmospheric paths can be considered the same), these 
problems can be solved by observing both sources simultane- 
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ously, or by rapidly switching between them, enabling their rel- 
ative phase to be tracked as the Earth rotates, and solving for 
the 2n phase ambiguities (Marcaide & Shapiro 1983). However, 
for very large source separations the relative phase between two 
sources cannot easily be used directly for relative astrometry. 

1.2. Group-delay astrometry 

Observations over a wide frequency band can solve the phase 
ambiguity problem by furnishing another variable; since (§> varies 
linearly with frequency for non-dispersive (geometric and tropo- 
spheric) signal paths, the phase slope across the band - the group 
delay t = (l/2n)d(f>/dv - provides an unambiguous measure of 
the path length difference, Lcos9/c, albeit with a reduced preci- 
sion which depends on the bandwidth spanned. A small position 
offset results in an astrometric delay change 

Ar flJfr = (l/2n)d(A<f> astr )/dv = LsinOAx/c. 

The dispersive ionospheric contribution produces a delay of 

T io „(v) = (h B -h A )v- 2 /c. 

By observing simultaneously in two widely separated frequency 
bands v\ , V2, the delay difference between the bands can be used 
to estimate, and remove, the ionospheric contributions to the 
group delay measurement at the higher frequency v 2 . 

At, „(v 2 ) = (t V2 - T Vl )/(l - (vi/v 2 r 2 ) . 
Group delay astrometry measurements of compact radio sources 
(mostly Active Galactic Nuclei, AGN) are regularly carried out 
by the International VLBI Service (IVS) in 24-hour observing 
runs, using a global network of telescopes observing at 2.3 and 
8.4 GHz, and bandwidths spanning a few hundred MHz. After 
ionospheric delay correction, the 8.4 GHz group delays are anal- 
ysed using a model for the tropospheric delay, producing fits 
for Earth rotation parameters, telescope coordinates and radio 
source positions. Such group delay positions of radio sources, 
and with quoted errors as small as 250 //as, are used to define the 
ICRF to an accuracy of ~30 //as. 

1.3. Astrometry with extended sources 

The visibility phase of an extended source is conventionally di- 
vided into 2 terms: (a) that due to the path length difference to 
some reference position within the source, x , and (b) a source 
structure term, (f> str , which can be determined by evaluating the 
relevant component of the Fourier transform of a map of the 
source brightness distribution centred on xo. Note that (f> str is 
generally a "baseline-dependent" quantity, which allows these 
2 terms to be separated by self-calibration procedures when 
imaging sources with a VLBI array (Cotton 1979; Cornwell 
& Wilkinson 1981; Schwab & Cotton 1983). All the sources 
used for defining the ICRF exhibit a dominant compact fea- 
ture - the "core" - which provides an obvious reference posi- 
tion within the structure. The difference of visibility phase be- 
tween two extended sources can be used to determine the pre- 
cise separation between their reference points by first imaging 
the sources, evaluating the structure phase terms and correct- 
ing their visibility phases (Marcaide & Shapiro 1983). Chariot 
(1990) has emphasised that this correction may also be nec- 
essary even when using group delays for astrometry, in those 
cases where <j> str varies rapidly with resolution, giving rise to 
a baseline-dependent structure delay (l/2n)d(f> str /dv across the 
observing band. In principle, source structure corrections should 
be made separately at both 2.3 and 8.4 GHz when ionospheric 
corrections are made to group delays (Petrov 2007). 



When sources are observed using a wide bandwidth, the 
change of source structure with frequency across the band needs 
to be taken into account. The technique of Multi-Frequency 
Synthesis (MFS) has been developed for astronomical imaging 
of extended sources, which uses wide (30%) bandwidths in or- 
der to fill in the range of resolutions covered by interferometer 
arrays (Conway et al. 1990; Sault & Wieringa 1994). The algo- 
rithms used for MFS are primarily concerned with producing a 
single-frequency image with high dynamic range, but they also 
estimate the spectral index distribution. Evaluating the source 
structure contribution to the wideband visibility function used 
for group-delay measurements in principle requires both of these 
components. 

2. Frequency-dependent core-shifts 

The sources used by the IVS for defining the ICRF are chosen to 
be largely point-like at both 2.3 and 8.4 GHz, making structural 
corrections apparently unimportant. However, due to opacity in 
the dominant, compact core component at the base of the rela- 
tivistic jet, the position of the peak of radio emission is expected 
to be frequency-dependent (Blandford & Konigl 1979; Konigl 
1981). Marcaide et al. (1985) discovered a shift of 700 //as be- 
tween the positions of the "cores" of the quasar 1038+528A 
measured at 2.3 and 8.4 GHz by comparing their separations 
from a feature in the nearby quasar 1038+528B. Using obser- 
vations at additional frequencies, they described the position 
change with an ad hoc law, kA^ 3 , with 0.7 < /3 < 2.0. In order 
to measure such "core-shifts" it is necessary to correctly reg- 
ister the images at different frequencies. This can be done by 
assuming that another feature of the structure has a frequency- 
independent position, by phase-reference imaging using another 
source with a frequency-independent position, by using the tech- 
nique of frequency phase-referencing (Middelberg et al. 2005; 
Rioja et al. 2005) or by comparing the separations of multiple 
gravitationally-lensed images of a single source at different fre- 
quencies (Porcas & Patnaik 1995; Mittal et al. 2006). 

Lobanov (1998) has suggested the application of such opac- 
ity effects to the study of the conditions within the central re- 
gions of AGNs. He measured a value ft = 1.04 for the core 
of the quasar 3C 345, close to the value of 1 expected for syn- 
chrotron self-absorption in the regime of equipartition between 
jet particle and magnetic field energy densities. Kovalev et al. 
(2008) have investigated the prevalence of core-shifts in a sam- 
ple of 29 sources for which a secondary (presumed achromatic) 
feature in the jet could be used as a position reference. They re- 
port a median core-shift of 440 /ias between 2.3 and 8.4 GHz, 
with a largest value of 1400 //as. Although the sources used in 
these studies necessarily had jet components in addition to the 
compact, "chromatic" core, there is no reason to suppose that 
such core-shifts are absent even in sources which show little evi- 
dence of other structure. Thus one must assume that the sources 
used to define the ICRF (at 8.4 GHz) in fact have frequency- 
dependent positions. The measured values of core-shifts are typ- 
ically smaller than the beamwidth of VLBI arrays used to image 
sources and this effect does not, therefore, easily lend itself to 
correction using the MFS techniques mentioned above. 

2. 1 . Effect of core-shifts on astrometry measurements 

Hitherto the effects of AGN core-shifts within the observing 
band on VLBI astrometric group delay measurements have not 
been widely considered. Here this effect is investigated. For sim- 
plicity the core is considered to be a point source at each fre- 
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Table 1. Position shifts in //as assuming various values of 
for a source with measured core-shift between v\ =2.3 GHz and 
V2=8.4GHz of 440 //as, the median value found by Kovalev et 
al. (2008) 



p 


Ax 2 .3 


Ax 8 . 4 


(1-j8)Ax 8 .4 


Ax- 


^■^■phase-delay 


0.2 


1928 


1488 


1190 


29 


326 


0.4 


1088 


648 


389 


21 


281 


0.6 


814 


374 


150 


14 


239 


0.8 


682 


242 


48 


7 


201 


1.0 


606 


166 








166 


1.2 


558 


118 


-24 


-7 


134 


1.4 


526 


86 


-34 


-14 


106 


1.6 


503 


63 


-38 


-21 


80 


1.8 


487 


47 


-38 


-29 


57 


2.0 


476 


36 


-36 


-36 


36 



Notes: columns 2 and 3 give the deduced shifts from the jet base x at 
2.3 and 8.4 GHz; column 4 the "reduced" group-delay shift at 8.4 GHz; 
column 5 the "ionospheric" shift; column 6 the total difference be- 
tween positions measured at 8.4 GHz using interferometer phase, and 
ionosphere-corrected group delays. 

quency, whose position along a fixed direction depends on fre- 
quency. Following Marcaide et al. (1985) the shift, Ax(A), with 
respect to a reference position, xo (the jet base), is parametrized 
with a power-law, Ax(A) = kA@ which can be re-written as 

Ax' (v) = Ax'j 

where Ax'j is the shift in the resolution direction at unit fre- 
quency. The astrometric phase term then becomes 

A<p astr — (27Tv/c)Lsin0 Ax^P . 
Relative astrometry using interferometer phase measurements 
refers to the frequency-shifted core position or the average posi- 
tion of the core within that band. The astrometric contribution to 
the group delay, however, becomes 

A Tastr = (l/2ji)d(A<f> astr )/dv = (1 -frLsMAx^-Plc. 
Note that the astrometric correction from a group delay measure- 
ment responds to a "reduced" core-shift of (1 -/3)Ax l v^. In the 
special case of ft = 1, the group delay measures no core-shift at 
all but gives a position corresponding to the jet base, xo. In gen- 
eral, phase and group-delay positions will differ at any frequency 
by an amount fiAx\v^ . 

The presence of core-shifts interacts with dual-frequency 
group-delay corrections for the ionospheric path. If not taken 
into account, the delay difference between measurements at two 
frequencies introduced by the core-shift, 

AT astr (v 2 ) - Ar astr (vi) = (1 -p)Lsin0 Ax\ {v 2 ~ p - vC p )/c 
will be incorporated as part of the ionospheric correction, 
At; „(v2), resulting in a further (but typically smaller) shift Ax, „ 
in the measured position at v 2 : 

^ion = (! -/?)A*;(V 2 ^ - Vi^)/(1 - (VIM)" 2 ) . 

Of course, for the case of ft = 1, there is no measured core-shift 
and hence the ionospheric correction is properly applied. 

3. Discussion 

Table 1 compares values of core-shifts and group-delay mea- 
surements for various assumed values of /3 for a source having a 
core-shift between 2.3 and 8.4 GHz of 440 //as, the median value 
found by Kovalev et al. (2008). For values of above 0.7 the 
typical shift from xo at 8.4 GHz is within the current quoted 
250/ias error for group-delay positions. For ft = 1 the typical 



shift is 166 //as. For a baseline length of 6000km (such as that 
between Westford, Ma, USA and Wettzell, Germany) and an ob- 
serving frequency of 8.4 GHz, Cliobe 1S 1-2 mas. The shift pro- 
duces an astrometric phase term A(f> astr = 49°. For a fixed offset 
of 166 //as the phase change across a 720 MHz band (as used for 
IVS observations) due to the changing resolution would amount 
to 4°. However, the small shift in core position across the band of 
14//as exactly cancels this change, resulting in a phase offset of 
49° at all frequencies and an astrometric group delay AT astr = 0. 
Note that, for f3 > 1, the group delay position is on the opposite 
side of xo from the jet direction. For the (probably physically un- 
likely) case of /5 = 2 the effect of the core-shift on group delays 
is indistinguishable from that of the ionosphere, and ionosphere- 
corrected delays again refer to xo, as in the case of /3 = 1 . 

3.1. The ICRF 

The ICRF is the current best realization of the Celestial 
Reference Frame, based on ionosphere-corrected, 8.4 GHz 
group-delay positions for several hundred radio sources. The 
sources are all compact and hopefully stable, resulting in consis- 
tent and repeatable positions when measured using the standard 
IVS observing bands. VLBI observations link the ICRF with the 
International Terrestial Reference Frame (ITRF) via group-delay 
positions of AGN cores in the sky and radio telescope positions 
on the Earth. Just as the mathematically well-defined, but some- 
times physically inaccessible, positions of telescopes (the inter- 
section of the azimuth and elevation axes) must be located with 
respect to the grid of geodetic monuments on the Earth, so too 
must "mathematical" group-delay positions of AGN cores be lo- 
cated with respect to radiating points in the sky at radio and 
other wavebands. The phenomenon of core-shifts complicates 
this procedure although, if /? = 1 for most AGN cores, the rela- 
tionship is less frequency-dependent than one would otherwise 
suppose, since all group delay positions at whatever frequency 
then refer to the same point. 

Jacobs & Sovers (2008) have recently established a Celestial 
Reference Frame based on VLBI group-delay measurements at 
8.4 and 32 GHz, where the 8.4 GHz measurements are used to 
make ionospheric delay corrections. Their frame is based on the 
positions of 318 compact radio sources. For sources in common 
with the ICRF, the weighted rms position differences are 241 //as 
in RA cos(dec) and 290 //as in declination. These differences are 
close to those one might expect from the measuring precisions 
alone. For the case of ft = 1, the group-delay positions at both 
frequencies measure the same position of the jet base, xo. Note 
that the typical core-shift at 32 GHz is only 44 //as, correspond- 
ing to the difference between phase and group-delay positions. 

3.2. Phase-reference observations 

VLBI phase-referencing is a technique whereby the visibility 
phase of one source (the calibrator) is subtracted from that of 
another (the target). It can be used either to remove unwanted 
instrumental and propagation phase terms from a weak target 
source visibility, and/or to measure the relative position between 
the target and calibrator. The technique has also been success- 
fully used for spacecraft navigation. Calibrator sources and their 
positions are often chosen from those used to define the ICRF, 
or from the VLB A Calibrator Survey (VCS; Petrov et al. 2008 
and references therein) for which similar precise ionosphere- 
corrected 8.4 GHz group-delay positions have been measured. 
As shown above, these positions do not, in general, refer to the 
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core position at 8.4 GHz. Although this does not detract at all 
from the ability of phase-referencing to detect weak sources, it 
does require a reinterpretation of the derived position of the tar- 
get, since this is measured with respect to the calibrator core 
position at the observing frequency, not its listed group-delay 
position. Knowledge of the calibrator core shift at the frequency 
of observation is hence necessary to locate the target in the ICRF. 

Relative positions derived from VLBI phase measurements 
are more precise than positions determined using group-delays. 
Fomalont (2006) has suggested that more accurate ICRF posi- 
tions might be obtained from 8.4 GHz relative phase measure- 
ments by first determining the relative positions of groups of 
sources within a radius of perhaps 20°, and then "stitching" 
together many such groups to cover the whole sky. This may 
indeed be possible but one should note that the relative posi- 
tions measured at 8.4 GHz between ICRF sources will not cor- 
respond to their relative positions measured using group-delays 
at 8.4 GHz. 

Marti-Vidal et al. (2008) have used wide-field phase- 
differenced observations at 15.4 GHz, and an atmospheric 
model, to determine positions for 13 AGN cores in the S5 po- 
lar cap sample, with quoted precisions ranging from 14-200 //as. 
They determined mean position corrections with respect to ICRF 
group-delay positions of 278 //as in RA and 170/ias in declina- 
tion. Note that, for/? = 1, the group-delay positions refer to the 
jet base, whereas the 15.4 GHz interferometer phase positions 
refer to the 15.4 GHz core position, typically ~93 //as away. 

3.3. VLBI2010 

Within the program VLBI2010 the IVS community is planning a 
major upgrade in instrumentation for global geodetic and astro- 
metric observing (Petrachenko 2006). Amongst other improve- 
ments it is planned to observe in 4 different frequency bands 
between 2 and 18 GHz, both to improve the group-delay res- 
olution, and also to attempt the transition from group-delay to 
using the full precision of interferometer phase measurements. 
Hobiger et al. (2009) have considered schemes for solving the 2n 
phase ambiguities which arise in interpolating the interferometer 
phase between the bands, taking into account the need to solve 
for the dispersive ionospheric path and core-shifts in the different 
bands. Their simulations, however, which use phase and group- 
delay measurements from the 4 bands, assume fixed core-shifts 
in each of the bands and do not take account of the shift within 
the bands. In fact, for the case ft = 1, the core-shift has no ef- 
fect on the group-delay measurements and adds just a constant 
phase offset at all 4 bands, which would make the interpolation 
problem in the presence of an unknown ionospheric path eas- 
ier to solve. However, for unknown values of /3 and Axi phase 
measurements in at least 3 frequency bands would be needed to 
determine these core-shift parameters. 



4. Conclusions 

It has been shown that the phenomenon of AGN frequency 
core-shifts, arising from opacity effects on the radio synchrotron 
emission at the base of relativistic jets, breaks the proportional 
relationship between interferometer phase and group-delay mea- 
surements in astrometry. The effect on the latter has been in- 
vestigated using a simple model of a point-source core whose 
position down the jet has a power-law dependence on wave- 
length. This is no doubt an over-simplification as both the ex- 
ponent p and the jet angle may change with position along the 



jet (Lobanov, 1998). Nevertheless, the results indicate that sig- 
nificant position differences arise at the sub-mas level between 
measurements using phase and group delays. These differences 
are comparable to the current precision of source positions de- 
rived from group delay measurements and used to establish the 
ICRF. However, they will assume a greater importance in the 
near future when improved radio-astronomical instrumentation, 
and the GAIA mission at optical wavelengths, produce positions 
with up to an order of magnitude better precision. 

The core-shift coefficient, Axi, and power-law exponent, 
/?, are governed by the precise physical conditions in the 
synchrotron-emitting jet. Blandford & Konigl (1979) have 
shown that, for synchrotron self-absorption in an equi-partition 
regime, Axi is proportional to (luminosity) 2 ^ 3 , whereas /3 can be 
expected to be ~1. If this is indeed common, then group-delay 
measurements at any frequency refer to a fixed point at the base 
of the jet, perhaps close to the centre of optical emission. This 
is more simple than one might expect. Note that the process of 
delay fringe-fitting, used for most astronomical VLBI data anal- 
ysis, determines delay (but not phase) residuals with respect to 
this position. An improved Celestial Reference Frame, based on 
interferometer phase measurements of core positions, would pre- 
sumably align well with the present ICRF, since the offsets be- 
tween the cores and the jet bases in the set of defining source will 
be randomly distributed in angle. The VLB 120 10 project holds 
promise for such an advance in the coming years. 

Acknowledgements. I thank Bill Petrachenko for discussions on the use of vis- 
ibility phase in future VLBI2010 observing scenarios, and Andrei Lobanov for 
comments on the text. 
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